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Abstract. 

A low energy theory of suspended carbon nanotube quantum dots in weak 
tunnelling coupling with metallic leads is presented. The focus is put on the dependence 
of the spectrum and the Franck-Condon factors on the geometry of the junction 
including several vibronic modes. The relative size and the relative position of the dot 
and its associated vibrons strongly influence the electromechanical properties of the 
system. A detailed analysis of the complete parameters space reveals different regimes: 
in the short vibron regime the tunnelling of an electron into the nanotube generates a 
plasmon-vibron excitation while in the long vibron regime polaron excitations dominate 
the scenario. The small, position dependent Franck-Condon couplings of the small 
vibron regime convert into uniform, large couplings in the long vibron regime. Selection 
rules for the excitations of the different plasmon-vibron modes via electronic tunnelling 
events are also derived. 
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1. Introduction 

The nanoelectromechanical systems (NEMS) are characterized by a peculiar interplay 
between electronic and mechanical degrees of freedom [1]. Suspended carbon nanotubes 
constitute a particularly interesting realization of NEMS due to their remarkable 
electronic and vibronic properties [21 |3l H]. NEMS can be realized, though, in a 
variety of different flavors, with single molecule junctions O EJ [TJ [H |9], suspended 
and laterally confined two dimensional electron gases [lO], silicon [11], or suspended 
graphene [12]. Interesting vibrational effects in electronic transport have been observed 
in several recent experiments on suspended single wall carbon nanotube (SWCNT) 
quantum dots [131 El ESj- These experimental works have triggered several attempts 
[T6l [TS] \T7\ [T8] to theoretically explain some characteristic features of the measured 
stability diagrams {i.e., of the differential conductance in a bias voltage-gate voltage 
colour map). In particular, the height of the conductance peaks associated with the 
vibronic resonances is in quantitative agreement [13] with predictions of a simple 
Franck-Condon model for a single electronic level coupled to a harmonic mode (the 
so called Anderson- Holstein model) [T9l l20l 1211 [22l |23] . Nevertheless, the size of the 
electron phonon couplings required to fit the experimental data has remained essentially 
impossible to achieve with a microscopic theory [161 without introducing large 
screening effects [ISj. Moreover, the experimental data present negative differential 
conductance features which go beyond the capability of the simple Anderson-Holstein 
model. Different extensions of this model [211 |25l ESI I2Z] including asymmetric coupling 
or multiple electronic levels have been proposed to explain NDCs. In a recent work [18], 
they have been attributed to a spatial dependent Franck-Condon factor, as it naturally 
occurs in a clamped nanotube, combined with the assumption of a vibron mode being 
mostly localized near one of the two dot ends. 

Convinced of the fundamental importance of the geometrical configuration of the 
junction on the transport characteristics of a suspended SWCNT, we improve and 
extend here the work presented in [18]. Specifically, we calculate the spectrum and 
the tunnelling matrix elements over the entire parameters space obtained varying the 
relative length and relative position of the vibron with respect to the quantum dot 
also including the effect of higher vibronic modes. The treatment of a wide parameters 
space is relevant since it allows for a unified picture of different results presented in 
the literature [13 [18]. Also the inclusion of the higher vibronic modes seemed to us 
a necessary extension for two reasons: there is no real energy separation between the 
different vibrational modes since the frequency of the nth vibronic mode is just and nth 
multiple of the fundamental frequency u of the lowest one; furthermore, the very same 
dispersion relation (linear with respect to the mode number n) implies the presence 
of several degenerate vibronic configurations for the system, a necessary condition for 
interference triggered NDC features in the stability diagrams of nanoj unctions in the 
single electron transistor set-up [23 1221 [30l El [32] • 

The spectrum is obtained via the exact diagonalization of the system Hamiltonian 
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with the mechanical degrees of freedom being coupled both to the total charge and the 
plasmons of the nanotube. As a result, at low energies the system is described by a set 
of displaced vibron-plasmon excitations and the tunnelling matrix elements reduce to 
the product of Franck-Condon factors, one for each vibron-plasmon mode. Importantly, 
the Franck-Condon couplings are different for the different modes and depend crucially 
on the geometry of the system. A detailed analysis of the complete parameters space 
reveals different regimes: in the short vibron regime the tunnelling of an electron into 
the nanotube generates a plasmon-vibron excitation while in the long vibron regime 
the polaron excitation dominates the scenario. The small, position dependent Franck- 
Condon couplings of the small vibron regime convert into uniform, large couplings in 
the plasmonic case. 
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Figure 1. Different realizations of a nanojunction with a suspended SWCNT. The 
lengths id of the quantum dot and of the vibrons are also indicated together with 
their position. The length and position of the vibrons is assumed to coincide with the 
suspended part of the tube. The length and position of the dot depend instead on 
many factors like for example the weak or strong hybridization of the SWCNT and 
the metallic leads and the presence of impurities or of side gates. The labelling of the 
different configurations is given according to the general one used in figure [2] 



The paper is organized as follows: In section 2 the model Hamiltonian of a 
suspended SWCNT coupled to several stretching modes is introduced. Particular 
emphasis is given to the dependence of the electron-vibron coupling on the geometry of 
the system, see figure [H A detailed analysis of the coupling constants as a function of 
the geometrical parameters is performed. A set of canonical transformation including a 
polaron transformation is employed in section 3 to obtain the spectrum of the SWCNT 
in the presence of electron-electron and electron-vibron interactions. Both analytical 
results for limiting cases and general numerical findings (see figure H]) on the entire 
parameters space are discussed. As known from the theory of Franck-Condon blockade in 
the simplest Anderson-Holstein model [201 122], the polaron transformation also crucially 
affects the tunnelling Hamiltonian describing the coupling to the source and drain leads. 
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In section 4 an analytical expression of the tunnelling matrix elements is derived. A 
detailed analysis of the associated Franck-Condon couplings is performed, revealing 
different regimes and selection rules for the tunnelling processes depending on the 
geometrical configurations. Conclusions are drawn in section 5. 

2. Low-energy Hamiltonian of suspended SWCNTs 

The low energy spectrum of finite size, interacting metallic SWCNTs has been discussed 
in [33] within a mean field approach and in [311 ES] within a bosonization framework 
going beyond the mean field results. Bosonization is also the natural approach to include 
the effects of the coupling to the longitudinal stretching modes [T71 [18]. Here, following 
[T7[ [T8| [M| . we derive and discuss the spectrum and many body states of suspended 
metallic SWCNTs. Particular emphasis will be given to the dependence of the electron- 
vibron coupling on the geometrical configuration of the system. An analytical expression 
of the electron- vibron coupling constants in terms of the relevant geometrical parameters 
is derived and plays a crucial role in the analysis of the spectrum and the matrix elements 
conducted in the following sections. 

We thus consider a Hamiltonian of the form 

-f^sys = Ho + Vcc + + Hcv, (1) 

where Hq is the noninteracting Hamiltonian of a finite size, metallic SWCNT, V^e 
describes the electron-electron interaction, is associated with the longitudinal 
stretching modes while Hev describes the electron-vibron coupling. 

2.1. Metallic nanotuhes at low energies 

Exemplarily we shall perform the quantitative analysis for armchair SWCNTs. The 
extension to arbitrary chiralities, though, does not change neither the essence of 
the calculations nor the main results presented here. In armchair SWCNTs at low 
energies and under periodic boundary conditions only the gapless energy subbands 
with linear dispersion touching at the Fermi points F = ±Koex (also denoted Dirac 
points), where ix is directed along the tube axis, are relevant [361 EZ]- Imposing open 
boundary conditions along the tube length Lj, the eigenfunctions of the noninteracting 
Hamiltonian Hq are the standing waves [311 ES] 

(r) = 4 E sgn(F)e-s'^(^)'^- Yl Ur^r>^ (0 , (2) 

where (ppp (r) describes fast oscillating Bloch waves on the graphene sublattice p = ± at 
the Fermi point F. The branch index r = ± accounts for left (+) and right (-) moving 
electrons. For armchair SWCNTs it holds /pr = if P = + and fpr = —^'^^P=—- 
The parameter 

TT 1 

K = — K + A), n^eZ, |A| < -, (3) 
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measures the wave number with respect to the Fermi points K^, while A has to be 
introduced if there is no integer n with Kq = j^. Including the spin degree of freedom 
cr the Hamiltonian Hq therefore reads 

Ho = hvp ^^Y^ i^cl^A'^K. (4) 

where vp ~ 8.1 x lO^m/s is the Fermi velocity. Thus the level spacing of the 
noninteracting system is given by eo = ^'^f^^ while = EqA denotes the energy 
mismatch between the r = ± branches. Moreover, the operator CraK annihilates an 
electron in the state {(prK)^). In turn, the electron field operator is expressed in terms 
of the wave functions iprn (0 as 

(r) = ^ (fru (r) CraK- (5) 

rK 

The electron-electron interaction assumes the standard form 

Ke = ^ XI / df'¥^{r)¥^,{f")U{f- f ')^a(r1, (6) 

cr,cr' 

where for the actual calculations we model U [f—f') by the so called Ohno potential [35] 



U{f-f') = Uo 



f/nelr — r " ^ ^ 



a 



-1/2 



(7) 



where a reasonable choice of the onsite energy is [35] Uq = 15 eV, the dielectric constant 
is e ^ 1.4 - 2.4 and a = 14.397 eVA. 

The Coulomb interaction causes Umklapp, backward and forward scattering 
processes among the electrons. Away from half filling it is reasonable to neglect Umklapp 
scattering. We also disregard backscattering processes, which is a valid approximation 
for nanotubes with not too small radii [35]. The forward scattering processes can be 
fully included within a Tomonaga Luttinger (TL) model for SWCNTs [TTJ |36] yielding 
the TL Hamiltonian: 

Ho + Kc ~ ^TL = i^N + (8) 

j 

where H^ describes the fermionic configuration of the nanotube and Hj represents the 
bosonic excitations with the index j = c+, s+, c— , s— labeling the four excitation sectors 
for total charge, total spin and relative (with respect to the two electronic subbands) 
charge and relative spin, respectively. The fermionic component of ([8]) can be casted 
into the form: 

^N = f Ef +^AiV.„+Ee^ (9) 

j 

where the particle number operators for the different charge and spin sectors are 
defined by Nc+ = Y^ra^ra, N^- = Em sgn(r)iVr^, Ns+ = Era sgn(o-)iV,.^ and 
Ns_ = 'Y^^sgn{ra)Nra, and the operator Nra counts the particles with spin cr and 
pseudospin r. The electron-electron interaction is parametrized, in the fermionic part of 
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the Hamiltonian, by E^, i.e. the charging energy of the SWCNT quantum dot. Finally, 
Hj describes the bosonic excitation of the sector j. In the long wavelength limit, it 
reads: 



where the sum runs over the mode number n. Due to the Coulomb interaction the factor 
Qj < 1 for the sector c+ while Qj = 1 for the other cases. For unscreened interaction 
^ 0.2 [351 [36]. 

2.2. The electron-vibron Hamiltonian 

The low-energy vibrational excitations of the nanotube can be described in terms of low- 
energy acoustic modes (3l[38l[39]. These modes are coupled to the electronic degrees of 
freedom either via a deformation potential (associated with local variations in area) or to 
bond lengths modifications. The latter coupling mechanism has a coupling constant one 
order of magnitude smaller than the one associated with the deformation potential [3] . 
Hence the twisting modes which involve pure shear and thus a modification of the 
bond length can be neglected. Likewise the bending and breathing modes, though 
coupled via the dominant deformation potential, do not play a significant role [31 [38]. In 
fact, the bending modes only couple quadratically to the electronic degrees of freedom, 
while the breathing modes lie too high in energy to be excited in low-bias transport 
experiments. Thus, in doubly clamped SWCNTs, the stretching modes only can be 
retained, in agreement with experimental conclusions jl3j . 

Following Ref. [17] the stretching mode Hamiltonian is expressed in a continuum model 

as 



where ( = 27rRM, with R the tube radius, M is the carbon mass per unit area and Vgt is 
the velocity of the longitudinal stretching mode. Moreover, and are the position 
of the center and the length of the vibron, respectively. Typical SWCNT parameters 
are: Vst = 2.4 x lO^m/s, M = 3.80 x lO'^kg/m^. 

Notice that the positions and the lengths of the dot and of the vibron do not 
necessarily coincide. The length of the vibron (L^) is readily estimated as the distance 
between the electrodes which clamp the nanotube and it is defined as the length of 
the free standing portion of the tube. Instead, the relation between the size of the 
quantum dot and the geometrical properties of the junction is much more complex. The 
best way to estimate is to extract it from transport measurements which give the 
mean level spacing and the charging energy of the system. The position of the center 
of the dot Xd can be taken free parameter. 

In figure [T] we sketch four possible physical realizations of different configurations. In 
the panels A and A' the dot lies inside the vibrating part of the tube. The confinement 
is obtained by a side gate (A) or by impurities located on the tube (A'), while the 




(10) 




(11) 
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rest of the tube is electrically absorbed into the leads due to the strong tube-lead 
hybridization (extended lead configuration). In panel C the dot coincides with the 
entire tube length due to the weak hybridization between the SWCNT and the metallic 
leads and fully contains the vibrating fraction of the tube. Finally, a somehow mixed 
scenario is illustrated in panel B. 

The electron-vibron coupling Hamiltonian reads: 



Hp 



dfp{f)V (f) 



(12) 



where p{r) = Xlo- ^3t(0^o-(0 is the electron-density and V{r) = gdxu{x) is the 
deformation potential for the stretching vibron mode. The coupling constant g is 
estimated to be [3] 5^ ~ 20 — 30eV. The displacement and momentum field operators 
read BOl 



ulx) 



P{x) 



^sin 



m>l 



E 

m>l 



sm 



X 



x^ + 



2 



(13) 



km (' 



with km = mvr/Lv the wave number. Here 0^(0^) the annihilation (creation) 
operators associated with the mth vibron mode obeying the commutation relation 



Using the above field operators, we obtain 
2 



m>l 



m I d^Clm ~\~ „ 



(14) 



with Em = mhVst'n'/L^ = mhw. Similarly f ll2p . integrated over radial and azimuthal 
coordinates, becomes 



Hcv = g 



m>l 



1/2 



(al^ + am) / dxpiD(a;)cos 
^dnv 



km I X Xv "t" 



(15) 



where the integral is calculated over the overlap of the dot and vibron region and the 
effective ID density operator pmix) reads, in its bosonized form [Mj : 

- . - ^ 2 



PlD[X) 



-.dx(l)c+{x). 



(16) 



Notice that the bosonic field (j)c+{x) can be expressed in terms of the bosonic creation 
and annihilation operators bl_^_ „ and 6c+,n as: 



0c+(a;) 



hg, 



c+ 



^sin 



n>l 



11 , Ld 

Kn \ X Xd -\ — 



Hi. 



(17) 



where kn = rnv/Ld and we have imposed open boundary conditions p{xd — Ld/2) = 
p{xd + Ld/2) = 0. It is useful, for the diagonalization procedure presented in the next 



Spectrum and Franck- Condon factors of interacting suspended SWCNTs 



8 



subsection, to introduce the dimensionless position and momentum operators for the 
nth plasmon mode {Xn, Pn} and the ones of the mth vibronic mode {xm,Pm} 



Xr, 



bc+,n + h\ 



c+,n 



V2 



-'c+,n 



Pv 



\V2 ' "™ iv^ ' 

which satisfy the canonical commutation relations [X„,P„/] = i(5„,„/ and [xmiPm'] = 
i6mm'- In terms of these operators the electron- vibron Hamiltonian can be written as: 



n,m>l 



(19) 



m>l 



where 



9\ 



(20) 



is the fundamental coupling constant; it acquires the value / = 88 fieV for a (10, 10) 
SWCNT with Ld = 1 fim and assuming g = 30 eV. 

The geometric part of the electron-vibron coupling is given by the dimensionless 
matrix 



-f^nm(A, S) 



dx < cos 



cos 



TTX (n + 



TTX n 





mvr 


j) 




m\ 


mvr 

+ - 


j) 



1 - A 
2 

1 - A 



(21) 



for the plasmon-vibron component and by the vector 



dx cos 



mvr 

— ix-6 



(22) 



(23) 



for the charge- vibron component. The integration limits 

Xmin = max[0, 5 + (1 - A)/2], 

a;i„ax = min[l,5+(l + A)/2] 

ensure that the integral extends only on the overlap regions of the dot and vibron. 
As one appreciates from (fT9l) -( !22l) . for a fixed the electron-vibron Hamiltonian is 
completely determined by the relative position of the centers of the dot and the vibron 
S = (xv — xa)/Lai and the ratio between the length of the vibron and of the dot 
A = L^l Lj. 

Importantly, ifev reveals that the electron-vibron interaction only involves the 
position operator Xm of the m-vibron mode, the position operator X„ of the n-th 
charged plasmon mode and the total electron number Nc+- Moreover, the important 
energy scales involved in the electron-vibron dynamics are the lowest vibron energy 
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huj, the lowest charged plasmon energy eo/gc+ and the fundamental coupling constant 
/: their values are O.ObOmeV, 8.293meV^ and 0.088meV , respectively, for a (10,10) 
SWCNT with Ld = = Ifim and the other parameters as the ones given in figure ID 
Excluding the extreme short vibron regime A < 1/100 and the strong screening, we can 
conclude that hw,I <^ ^o/9c+ thus implying a clear separation of the vibron and plasmon 
energy scales. Albeit these two degrees of freedom are consequently characterized, for 
an isolated system, by an essentially independent dynamics, the tunnelling event can 
be substantially influenced by the mechanical motion of the nanotube under certain 
geometrical conditions as will be discussed later. 

2.3. Plasmon-vibron and charge-vibron couplings 

The energy spectrum and the Franck-Condon couplings strongly depend on the geometry 
of the junction via the coupling constants Knm and L^. The detailed analysis of 
these coupling constants is thus the natural starting point to understand the presence 
of geometrical dependent trends and selection rules in the tunnelling processes of a 
suspended SWCNT. 

The geometrical parameters space {A, 6} is divided into four regions by the different 
conditions imposed by the integration limits Xmin and Xmax explicitly given in ( l23l) . In 
figure [2] we define these regions and give a schematic representation of the corresponding 
geometrical configuration. 




1 



Figure 2. Parameters space of the geometrical configurations of tlie electromechanical 
nanojunction. The relevant dimensionless parameters are the length ratio A = L^/L^ 
and the relative position of the centers d = (xv — Xd)/Ld. Four qualitatively different 
regions are identified in the parameters space and schematically shown on the right. 
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The function K„,m has the following explicit form in the four regions: 



2 m- 



2 



-l)"sin 



1 - 2(5 + A 



niTT- 



2\n 



l)"m sin I mn 



2A 

1 - 2(5 + A 



26 -X 



sm mTT 



2A 



7r(A^?7,2 — m^) 
2 



;-l)'"sin Xrni 



2A 

A + 25 + 1 



-l^^Ansin Xmr 



2A 

A + 25 + 1 
2A 



An sin ( Xmr 



sin ( Xnn 



1 + 26-X 



+ m sm miT 



2X 

X-26-1 
2A 

A + 25 - 1 
^A 



(25) 



7r(A^n^ 

Some symmetry relations for the function Knm can be readily obtained from ( !24l) : 

i^„^(A,5) = l/Ai^„„(l/A,-5/A), 

i^nm(A,5) = (-l)"+"K„„(A,-5). 

The first equation in ( l25i) quantifies the connection between the behaviour of Knm{X, S) 
at small and large values of A: the roles of the vibron and of the plasmon are simply 
exchanged in the plasmon-vibron Hamiltonian if we invert the ratio of their lengths. 
The second equation in fl25l) states instead that if we invert the relative position of the 
vibron and the dot, the plasmon-vibron Hamiltonian acquires at most a minus sign, 
depending on the parity of the vibronic and plasmonic modes. The case considered in 
[17] corresponds to the point A = 1, 5 = of the parameters space where the following 
limit holds: 



lim KnmiX, 0) = 5r, 



(26) 



and each vibronic mode is only coupled to the plasmonic mode of the same order 
n = m. In all other regions of the parameters space the coupling is not diagonal 
and the dynamics of each vibronic mode is influenced by all plasmonic modes and vice 
versa, making the system quite intricate. Nevertheless, from a detailed analysis of the 
K function, one can estimate which modes are more relevant in the low energy limit. 

The function K has an upper bound K < 2, as it can be easily proven from its 
definition ( 1211) by considering that the distance between the integration limits is at 
maximum A. Thus, Knm does not diverge for A — )■ m/n as one could expect from a first 
sight. Instead, its maximum can be estimated by calculating the limit A — )■ m/n. One 
obtains: 



lini Kmn{X,6) 




for n < m, 
for n > m. 



(27) 



where the first and the second case are calculated in region A and C, respectively (see 



figure [2]). The absolute value l-ft'nml of the coupling exhibits |m — + 1 local maxima 
as a function of the relative displacement 6, separated by nodes in which the m-th 
vibronic mode is decoupled from the n-th plasmonic one. Notice that in the limit 
A — m/n the wavelength of the m-th vibronic mode coincides with the one of the n-th 
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plasmonic one, thus giving a physical interpretation to the resonance. One concludes 
that each geometrical configuration optimizes the coupling between specific plasmonic 
and vibronic modes. Moreover, the coupling between low vibronic modes and higher 
plasmonic ones is reached for short vibrons and is more efficient than the coupling of a 
low plasmonic mode with higher vibronic ones obtained, instead, for large vibrons. 

Another interesting regime can be reached in the small vibron region when the 
center of the vibron lies in the vicinity of the border of the dot. Let us consider for this 
reason the function Knm in the region B and with A ^ m/n. The following relation 
holds: 



Knm (\,l + ax] = —{-IT sin 



1 

mvr I a 

2 



, ____ ..... , , (28) 

where |a| < 1/2. The absolute value l-ft'nml of the coupling exhibits m local maxima 
as a function of a in the region B which are independent of the plasmonic mode 
n. This specific configuration has been chosen in [18] to describe a system in which 
the renormalization of the lowest vibronic mode due to the coupling to the plasmons 
produces a strongly inhomogeneous Franck Condon coupling in the tunnelling matrix 
elements to the carbon nanotube. In order to illustrate the arguments presented so far, 
we plot in figure ([3]) the plasmon-vibron couplings Ki^ and K^i as a function of A and 
6. Clearly visible are the maxima of the coupling close to A = m/n and the fan shape 
structure of the coupling close to the points {0, ±1/2} more visible in the case Ki^ due 
to the conditions given above. 

The second line of the electron-vibron Hamiltonian ( lT9l) describes the charge- vibron 
interaction and is proportional to the function Lm defined in ( l22l) . Also the coupling 
Lm, is defined on the parameters space {A, 6} by different functions in the four different 
regions already introduced for Knm- 

l-26 + X\ f 1 + 2(5 -A 

sm I mvr ^ + sm mvr 





1 


mvr 




1 


mvr 




0, 




1 
mvr 



2A / V 2A 

1-25 + A^ 



sm mvr 



2A J (29) 
1 + 25- 



sm mvr ; 

2A 

A symmetry relation can also be derived for this coupling, namely: 

LM,S) = i-l)"'LM.-S)- (30) 

The function Lm vanishes identically in the region C thus implying no charge-vibron 
coupling for systems in which the vibron is entirely contained inside the dot. The finite 
local coupling is in fact averaged away by the sinusoidal shape of the vibron. The form 
of Lm in the region B is readily understood with the parametrization: 



^ 1 x\ 1 

A, ± — h aA = sm 



2 / mvr 



1 

mvr I a 

2 



(31) 
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Figure 3. The plasmon-vibron Knm and the charge-vibron L„i couphng constants are 
plotted in the geometrical parameters space. Top row: On the left (right) the coupling 
between the first (fifth) plasmonic and the fifth (first) vibronic modes. Bottom row: 
Examples of the charge-vibron coupling Lm are given for the first (left) and the second 
(right) vibronic mode. Black solid lines indicate in all figures the borders of the regions 
A, B, C, D indicated in figure ([2]) and explained in the text. 



with \a\ < 1/2. As shown in figure [3l \Lm\ has indeed in the small vibron limit a 
fan shape with m maxima of magnitude l/(m7r) separated by m — 1 nodes. Thus, 
the charge-vibron coupling decreases for higher vibron modes and is also very sensitive 
to the geometry of the system. The geometry of the system even introduces selection 
rules: for example for a system with A = 0.1 and 6 = 1/2, (a = 0) only odd modes 
(m = 2a — 1, a G N^) exhibit a charge-vibron coupling. 

The maximum charge-vibron coupling for the mode m is reached when \Lm\ = 
2/(m7r) and is obtained only for vibrons larger than the dot (A region, A > 1) and 
centered with respect to it {6 = 0). Only even vibronic modes couple to the charge 
if 5 = and \Lm\ exhibits m/2 maxima in the positions A = m/{l + 2r) where 
< r < m/2 - 1, r e N. 

In conclusion, the electron- vibron coupling is very sensitive to the geometry of the 
junction, both in its plasmon-vibron and charge-vibron components. In general, for 
more symmetric systems {6 ~ 0) the plasmon-vibron component dominates the short 
vibron limit (A < 1, region C) while the large vibron limit (A > 1, region A) is dominated 
by the charge-vibron interaction. Only for a strongly asymmetric system {6 ~ ±1/2) 
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in the short vibron hmit (A ^ m/n) the two components can have the same strength. 
Moreover, in general, the strength of the couphng decreases with the vibron mode. 

Yet, the position of the nodes of the functions Knm and Lm depends on the vibron 
and plasmon mode numbers n and m, generating selection rules that depend sensitively 
on the geometry of the system. 



3. Diagonalization and spectrum 



Because the electron-vibron coupling only involves the total charge operator Nc+ and the 
plasmon excitations, the part of the system Hamiltonian which is still to be diagonalized 
is: 

h^l / 'n\ hoj 



'mLmV2Nc+Xr, 



n>l 

n,m>l 



m>l 



(32) 



where we have introduced the frequencies Q = tivy / {gc+Ld) and u = nVst/L^- The 
exact diagonalization of the Hamiltonian ( 132|) can be achieved in two steps: i) a set 
of canonical transformations eliminates the plasmon- vibron component; ii) a polaron 
transformation eliminates the charge- vibron component. The final result is a collection 
of shifted plasmon-vibron oscillators. 

The first step in the diagonalization is better understood by setting the plasmon- 
vibron part of the Hamiltonian into a quadratic form: 



sys 



X 
P 













\ 
















X 








Hpp 







P 


V 










/ 


VP/ 



(33) 



V p / 

where the components of matrix Hm defining the quadratic form are given by: 

(^pp)nn' = nhn/2 6nn', {H^^)mm' = rnhw/2 6mm' and (i^pv)„^ = {H^p)mn = 

I-^ygc+y/nmRnm- Moreover, we have introduced the vector of operators X = 
[Xi,X2,...]^ and analogously for x, P, and p. Finally we have defined the charge- 
vibron Hamiltonian Hcv The quadratic form in ( l33l) is simplified via the following set 
of canonical transformations: the first is the contraction 



m 
I 



1 / V mhjjj Xr, 

Vmhwp^rn, 



(34) 



that transforms the momentum block of Hm into the matrix 1/2. Notice that the 
commutation relations between position and momentum operators are conserved for 
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each mode: [X!^,P^,] = i6nn' and = i6nn'- Afterwards we perform the rotation 



x' 

+171 "^m' 



n=l m=l 



n=l m=l /nr\ 



= 

n=l m=l 

that diagonahzes the position block of Hu written in the primed variables. We have 
also introduced the total number of vibron (plasmon) modes A'v (Np) This can be done 
without loss of generality due to the presence of physical cut-off's both for the plasmonic 
and vibronic mode numbers. This transformation is physically the most important since 
it generates the position and momentum operators C,i and which identify A'p + A^v 
mixed plasmon-vibron modes. The matrix defining the quadratic form reads, in this 
mixed basis: 



M 



A 








1/2 



(36) 



where A is a diagonal matrix whose diagonal element A; defines the energy of the 
plasmon-vibron mode hui = \/2Ai. This relation becomes clear after the last canonical 
transformation, the expansion 

^ !^ (37) 

that brings the system Hamiltonian into the form: 

^:ys = ET^(^"' + ^') + ^- (38) 
I 

The effect of the transformations 0341) . (15^ and (1371) on the charge- vibron Hamiltonian 
Hcv is readily obtained: 

^cv = lV2y2 . I^UNp+m, I ilNc+. (39) 

The presence of Hcv requires a second step in the diagonalization procedure: the 
polaron transformation i^'sys = e "^-f^sys^^"^ where 

S = iV2^^mLm^l^UNp+m,iTTiNc+ (40) 

Im 



yielding 




Ks = E ^te' + -^') - E ^ I j:LmJ^U,^+mJ ] Nl- (41) 
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Thus, the low energy spectrum of the suspended SWCNT reads 

^7V,m = ^A? + XI + ^) + X ^^Omn,j (42) 

where N = [Nc+, Nc-, Ns+, Ns-] is the vector defining the electronic configuration and 
Ej^ the associated energy as can be computed from and fHTl) . The vector m, 
instead, contains the occupation numbers mi of the plasmon-vibron modes and the 
ones {mnj, j 7^ c+) of the other relative charge and spin bosonic modes. 

The diagonalization procedure presented here reproduces known results in some 
limiting cases. In the totally symmetric case (5 = 0, A = 1) where length and center 
of the dot and vibrons coincide, only the coupling between plasmons and vibrons with 
the same number of modes is allowed {Knm = ^nm)- One obtains that the matrix to be 
diagonalized by the rotation (!35|) is: 



n' I X- I (43) 



yielding the spectrum [T7j: 

hioi = \/2Ai = nh\ 



n^ + u^^ Uil-u^Y^Ag^^ (44) 



where / = {n, a} and a = ±. For this symmetric configuration there is also no polaron 
shift since the charge- vibron coupling vanishes identically {Lm = 0). Also the case 
considered in [18] of a single vibron mode is reproduced by our general theory. Under 
the only assumption that u <^ fl one obtains: 



Ui = u 



00 



n=l 



which is always real for the parameters considered in the present paper. Moreover the 
case of short asymmetric vibrons (A <^ 1,5 = 1/2) is particularly interesting since by 
means of ( |28l) one obtains also that the lowest plasmons (n ^ 1/A) equally contribute 
to soften the frequency of the lowest vibron mode. 

In the generic case, though, only a numerical evaluation of the spectrum is viable. 
In figure m we present the relative frequency shift {i.e. {um — mu)/mu) for the first 
(left panel) and the fifth (right panel) plasmon-vibron mode. The calculation is carried 
out for a (10, 10) armchair nanotube of = l/im. The coupling of the vibrons to the 
plasmons softens the vibronic modes, yielding a negative shift for every configuration. 
The renormalization is stronger and almost constant in the region C where the coupling 
between the low vibronic modes to the plasmonic ones is larger. An estimate of the 
maximum renormalization can be obtained by its direct calculation in the symmetric 
point (A = 1, 5 = 0): 
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where we made the expansion of fl44l) in powers oi u/Q and 1/ {HQ). Interestingly, as far 
as the bare vibron frequency u and the fundamental electron-vibron coupling constant 
I/h are both much smaller than the bare plasmon frequency Q, the relative frequency 
normalization, if present, is independent of the mode number m. It is also clear that, 
in absence of strong screening, {gc+ ~ 0.2) the relative normalization is very moderate 
and does not exceed the 1 percent, independently of the geometry of the junction. 




Figure 4. Relative normalization of the first (left) and fifth (right) vibron mode due to 
the plasmon- vibron coupling. Notice that the relative normalization is always negative, 
it reaches its maximum in the point A = f , (5 = and is essentially constant in the entire 
region C (defined in figure [3]). The parameters used are Ld = 1 R ~ 6.68 A, v-p = 
8 X lO^m/s, v^t = 2.4 x W^m/s, M = 3.8 x IQ-^kg/m^, g = 30 eV, gc+ = 0.2. 



4. Tunnelling amplitudes and Pranck-Condon couplings 

So far we studied the isolated nanotube. Our interest, though, is the transport of 
electrons through a SWCNT in tunnelling coupling with possibly extended source and 
drain leads (see figured]). The tunnelling Hamiltonian Ht is given by: 

(47) 

a=s,d cr 

where see (jS]), and (r) = J2q^*qi'^ '^'^qaa electron creation operators in the 
SWCNT and in the lead a, respectively, and Tq {"f) describes the transparency of the 
tunnelling contact a. 

The spatial dependence of the transparency Tq (r) depends on the specific 
geometrical configuration of the junction and on the tube- lead hybridization. For the 
sake of simplicity we refer again to the configurations introduced in figure [H In both 
the A and A' cases we expect (r) to be strongly localized at the interface between the 
extended lead and the dot, while in the case C the tunnelling region extends over the 
entire fraction of the tube which is covered by the leads. For the case B an intermediate 
situation is envisaged with an extended tunnelling region (weak hybridization) at the 
source and a localized one at the drain (strong hybridization). 

In the weak tunnelling limit the dynamic of the system can be described as a series 
of sequential tunnelling events connecting different many-body eigenstates of the system. 



EE 



dr 



(r) (r) (r) + h.c. 
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For this reason a central role is played in the theory by the spectrum that we calculated 
in the previous section and by the tunnelling amplitudes between the corresponding 
many-body energy eigenstates which is the focus of the present one. 

Following [35] the 3D electron annihilator in the quantum dot ^^(r) can be 
expressed in terms of the slow varying ID operators iprFa{x) which, in their bosonized 
form, read: 



(4^ 



where 7)^0 is the Klein factor which reduces by one the occupation of the branch ro", 
KrFu{.x) is the operator 



KrFaix) 



1 i^sgn{F){rNrcT+A)i 



(49) 



V2L^ 

which essentially adds a phase proportional to the occupation number of the branch ra, 
and (j)rFa{x) is the bosonic field associated with the bosonic excitation of the SWCNT. 
It is useful to express iprFa{x) in terms of the position and momentum operators of 
the plasmonic modes and P„. After a lengthy but straightforward calculation one 
obtains: 



1prFa{x) OC Vra^rFaix) Y[e 

n>l 

where we have introduced the functions 

nn 



+iP„(x)X„-iX„(x)P„ 



(50) 



Pnix) 



COS 



ngc+ 



X 



n 



sgn(Fr) sin 



(51) 



and the proportionality in (l50l) is due the frozen c— , s+ and s— branches. They only 
contribute in fact with an overall constant to the tunnelling matrix elements between 
the low energy eigenstates. 

An explicit representation of these low energy eigenstates is readily obtained from 
(|9]) and ( 14T|) . Due to the already mentioned energy scale separation between on one 
side the vibronic and on the other side the plasmonic and electronic excitations, we can 
limit ourselves, without loss of generality, to the case m„^j = 0, j 7^ c+ and obtain 



where 



\N,rn)o = ll 



|iV,0)o. 



(52) 



(53) 



with = [Nc+, Nc-, Ng^, Ns_] being the vector defining the electronic configuration 
and m representing here the occupation numbers of only the lowest vibron-plasmon 
modes (with an excitation energy lower that Eq). The low energy eigenstates of a 
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metallic suspended SWCNT are, thus, polaron shifted plasmon-vibron excitations over 
its electronic ground state |A^, 0)o. We are now ready to evaluate the matrix element: 



{N,m\^rFa{x)\N',rn') = o{N,rh\e-^rFa{x)e+^\N' ,m')o. 
Since the operator S defined in ( HOl) commutes with KrFa- 



(54) 



(55) 



where the proportionality accounts for the constant terms deriving by the application 
of the Baker-Hausdorff theorem and we defined the functions: 

V2I 



^i{x) = > ^ \ — mLmUN^ 

£i 1 V £i 

m=l 



+m, I 



+E 



2e, 



Ur,.l COS 



(56) 



E 

n=l 



2gc+hQ 



Uni sin 



By means of ( 155|) it is now clear that the tunnelling matrix element factorizes into an 
electronic component and a product of Franck-Condon factors, one for each plasmon- 
vibron mode: 



{N,rh\iJrFa{x)\N',rn!) oc {N\firaKrFa\N')l[F{mi,m[, \i) 



where 



A, 



V2 



(57) 



(5J 



is the effective coupling between the charge and the plasmon-vibron mode and 
F{m,m', A) = \e{m' - m)A'"'^'" + e{m - m')(-A*)™-" 



my 



m 



(59) 



is the explicit expression of the Franck-Condon factor. The equations ( !56l) -( l59l) together 
with 0351) for the definition of the transformation U represent the main analytical result 
of this paper. They are a very general expression of the tunnelling matrix elements 
between the low energy eigenstates of a suspended SWCNT in presence of multiple 
plasmon and vibron modes. Special limits of these formulas are already present in 
the literature [T71 HB] . Particularly interesting to our point of view is the contribution 
of the geometrical configuration of the junction, which determines selection rules in 
the tunnelling processes and in turn the magnitude of the dimensionless electron- vibron 
Franck-Condon couplings A;. In figure [5] we present \Xi\ for the first (left) and the second 
(right) plasmon-vibron modes. The values in the figure correspond to a tunnelling matrix 
element calculated at the beginning of the tube [x = Xd — -f/d/S). By a comparison with 
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Figure 5. Franck-Condon couplings |A;| as a function of the length ratio A = L^/Ld 
and relative center position 5 = [x^ — Xd)/Ld. The coupling for the first and 
second vibron-plasmon modes are shown in the left and right panel, respectively. The 
parameters are the same as those reported in figure H) The couplings are calculated 
for a tunnelling event at the beginning of the dot. 



the corresponding charge-vibron coupling constant in the lower panels of figure ([3]) 
one can argue that \Xm\ oc \Lfn\. This observation is essentially correct, at least in the 
A, B and D regions of the parameters space where the energy renormalization of the 
vibronic modes is negligible and the same holds for the mixing introduced in equation 
(!35|) between the vibronic and plasmonic modes. Consequently, we expect that |A;| does 
not depend on the tunnelling point, at least in the long vibron region (A > 1) for any 
geometrical configuration. This result is illustrated in figure [6] where the Franck-Condon 
couplings for the first and second plasmon-vibron modes are plotted as a function of 
the dimensionless tunnelling point = x/La) and relative position of the vibron (5) for 
the configuration L^/L^ = 2. Interestingly, the selection rules derived in the previous 
section for Lm directly apply to the Franck-Condon couplings in the long vibron regime: 
for example for a symmetric junction {6 ~ 0) only even modes can be excited by a 
tunnelling event while the odd ones will remain in their ground state. Finally, it is also 
notable that max(|A;|) 1 in the long vibron regime even in absence of strong screening 
{9C+-0.2). 

A different result characterizes the short vibron limit (A < 1). In the C region 
the charge-vibron coupling vanishes identically due to symmetry considerations. Even 
if small, the vibron-plasmon mixing becomes there the dominant effect. In figure [7] 
we present the Franck-Condon coupling of the lowest vibron-plasmon mode for the 
configuration A = L^/L^ = 0.1. In particular, in the upper left panel we show |Ai| and 
in the remaining panels its components: i. e. in the upper right panel the charge-vibron 
component - the first line of C,i in fl56|l -, in the lower left panel the plasmon-vibron 
component of and vr; in the lower right panel. The dashed white lines represent in 
all panels the borders of the C region i.e. the region in which the vibron is completely 
inside the dot. Outside the C region the charge-vibron coupling is stronger and Ai 
does not depend on the tunnelling point. Inside the C region, instead, the dominant 
contribution is given (in the lower left panel) by the plasmon-vibron component of 
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Figure 6. Franck- Condon couplings of the first (left) and second (right) vibron- 
plasmon mode in the long vibron regime (A = 2) plotted against the dimensionless 
tunnelling point ^ = x/Ld and the relative position S of the dot and vibron centers. 
Dashed lines indicate the borders between the B (top), A (center), and D (bottom) 
regions of the parameters space (see figure [2]). 

The latter follows the position of the vibron and mimics its shape. The last observation 
is also confirmed by the lower left panel of figure |8] where the corresponding component 
of the Franck-Condon coupling for the fifth mode is plotted. Finally, for higher modes in 
the short vibron limit, a position dependent Franck-Condon coupling is still appreciable 
also in the B and D regions (see figure [HI). 

The relevance of these results for the tunnelling Hamiltonian and the associated 
tunnelling rates between the many-body eigenstates depends by their interplay with 
the spatially dependent transparency T(r) introduced in the beginning of this section. 
In fact we would expect to detect a position dependent Franck-Condon factor in the 
tunnelling rates only for the cases illustrated in figures lA, lA' and IB where the vibron 
extends also beyond the dot region but not for the case in figure IC. This observation, 
together with the results presented in figures EHHl allows to conclude that the position 
dependent rates can be observed, among the configurations considered in this paper, 
only in the asymmetric short vibron one (A ^ 1, 5 ~ ±1/2), i.e. a configuration 
of type B (or D), also in agreement with the results presented in [TBJ. In a recent 
publication [41j an alternative set-up has been proposed for the direct visualization 
of the position dependent Franck-Condon couplings in which one of the two metallic 
electrodes is substituted by the tip of a scanning tunnelling microscope. 

In absence of electron- vibron coupling the frequency of the nth stretching mode is 
an nth multiple of the frequency u of the fundamental mode. Hence, naturally, there 
are several energetic degenerate vibronic configurations (involving two or more modes) 
which may contribute to transport at finite bias. As we just proved, for realistic values of 
the parameters the softening of the stretching modes introduced by the electron- vibron 
coupling does not really lift these degeneracies. This fact has profound implications for 
the transport properties of the system. Interference effects have been in fact predicted 
even for systems in the Coulomb blockade regime [28l [311 1321 112] in presence of quasi- 
degenerate states. 
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Figure 7. Position dependence of the Franck-Condon coupling of the first vibron- 
plasmon mode in the short vibron regime (A = 0.1). In the upper left panel the full 
coupling I All is plotted while in the remaining panels its different components: i.e. 
in the upper right panel the charge- vibron component - the first line of ^/ in (|56p -, 
in the lower left panel the plasmon- vibron component of ^/ , and tt/ in the lower right 
panel. The dashed white lines represent in all panels the borders of the C region i.e. 
the region in which the vibron is completely inside the dot. 

Technically, this degeneracy determines the method of choice for the description 
of the dynamics of the system. At low biases, such that only the lowest vibronic 
mode is excited, a description of the dynamics only in terms of rate equations involving 
occupation probabilities of the many-body states of the quantum dots is appropriate. 
Yet, at higher bias, when several vibron modes are excited, a generalized master equation 
(GME) coupling diagonal (populations) and off-diagonal (coherences) elements of the 
quantum dot reduced density matrix should be used (see e.g. ([281 1221 SSI HS], HH HH 
SHI mi Ha So]). 

The sensitive dependence of the tunnelling matrix elements on the mode number 
for a given geometry of the system also suggests the existence of symmetrically coupled 
slow channels as the ones described in [27] and consequently of similar NDC effects in 
the stability diagrams of a suspended SWCNT junction. 

5. Conclusions 

In this paper we analyzed the spectrum and the effective Franck-Condon couplings of a 
suspended SWCNT quantum dot including many vibronic modes as well as different dot- 
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Figure 8. Position dependence of the Franck-Condon coupling of the fifth vibron- 
plasmon mode in the short vibron regime (A = 0.1). In the upper left panel the 
full coupling I As I is plotted. As in figure [7] the other panels represent its different 
components. 



vibron geometrical configurations. We described the long-wavelength acoustic-vibrons 
within an elastic continuum model and the electron-vibron interaction in terms of a 
deformation potential. In particular, we studied how the electromechanical properties 
depend on the relative size A and position S of the vibron with respect to the dot. 

Specifically, within the framework of a Tomonaga-Luttinger liquid description of 
the SWCNT, we derived an effective low energy Hamiltonian where the electron-vibron 
coupling is separated into a plasmon- vibron and a charge- vibron component proportional 
to different coupling constants {Knm and Lm, respectively). 

The system was diagonalized via a series of canonical transformations with an 
intuitive geometrical interpretation which reduce the low energy description of the 
suspended SWCNT to a set of displaced plasmon-vibron excitations. Consequently, 
the tunnelling matrix elements between the many-body eigenstates of the system are 
the product of Franck-Condon factors, one for each plasmon-vibron mode, of which we 
gave an analytical expression. 

The analysis of the coupling constants Knm and Lm and of the Franck-Condon 
couplings A/ on the entire geometrical parameters space allowed us to identify different 
regimes. 

In the short symmetric vibron regime (A < 1, 5 ~ 0) the charge- vibron component 
vanishes and the Franck-Condon couplings are extremely small {\Xm\ ~ 10~^) due to 
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the energy scale separation between the plasmonic and vibronic modes {Q/u ^ 1 and 
hQ/I ^ 1) that hinder the plasmon-vibron mixing. The Franck-Condon couphng is 
position dependent and is located around the position of the vibron. 

In the long vibron regime (A ^ 1) the charge- vibron coupling dominates 
the scenario giving substantially larger Franck-Condon couplings (|Am| ~ 1) and 
independent of the position as in the simple Anderson-Holstein model. The Franck- 
Condon couplings are strongly dependent on the relative position of the vibron and the 
dot, leading to selections rules: for example only even vibron-plasmon modes can be 
excited by electron tunnelling in a symmetric {6 = 0) long vibron junction (see figure 

In the asymmetric short vibron regime (A < 1, 6 ~ =tl/2) the charge- vibron 
and plasmon-vibron contribution are of the same order and correspondingly one 
can distinguish (at least in the higher modes, see figure [H]) the position dependent 
contribution due the plasmon-vibron mixing superimposed to the uniform polaron 
shift typical of the charge-vibron component of the coupling. In absence of screening 
(dimensionless electron-electron interaction strength g^^ = 0.2), though, the absolute 
value of the Franck-Condon coupling remains negligibly small compared to the one 
estimated from the experiments [131 [131 [ISj. Reasonable values have been obtained 
in this regime by [18] by assuming a very strong screening {g^^ ~ 1) that essentially 
removes the energy scale separation between the plasmon modes and the much shorter 
vibron mode. 

Finally, for reasonable values of the nanotube parameters the spectrum of the 
nanotube is only slightly modified by the electron-vibron coupling thus preserving the 
high degeneracy of the different vibronic configurations. This, in combination with the 
sensitive dependence of the tunnelling matrix elements on the mode number and on 
the geometry of the system, also suggests the existence of symmetrically coupled slow 
channels as the ones described in [27| and consequently of similar NDC effects in the 
stability diagrams of a suspended nanotube junction. 
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